function [f, gamma, eta] = rcnl_gmmobj_supply( ...
    kappa,   ...
    cdindex, ...
    cdid,    ...
    firmid,  ...
    prodid,  ...
    yearid,  ...
    bmc,     ...
    s_jt,    ...
    p_jt,    ...
    w,       ...
    derMat,  ...
    IVS,     ...
    W,       ...
    fes      ...
)

[eta, gamma] = rcnl_eta( ...
    kappa,   ...
    cdindex, ...
    cdid,    ...
    firmid,  ...
    prodid,  ...
    yearid,  ...
    bmc,     ...
    s_jt,    ...
    p_jt,    ...
    w,       ...
    derMat,  ...
    fes      ...
);

sMoments = eta' * IVS;
f        = sMoments * W * sMoments';
